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Abstract 

We consider Vlasov-type scaling for Markov evolution of birth- 
and-death type in continuum, which is based on a proper scaling of 
corresponding Markov generators and has an algorithmic realization in 
terms of related hierarchical chains of correlation functions equations. 
The existence of rescaled and limiting evolutions of correlation func- 
tions as well as convergence to the limiting evolution are shown. The 
obtained results enable to derive a non-linear Vlasov-type equation 
for the density of the limiting system. 
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1 Introduction 

The Vlasov equation is a famous example of a kinetic equation which de- 
scribes the dynamical behavior of a many-body system. In physics, it char- 
acterizes the Hamiltonian motion of an infinite particle system influenced by 
weak long-range forces in the mean field scaling limit. The detailed exposition 
of the Vlasov scaling for the Hamiltonian dynamics was given by W.Braun 
and K.Hepp [3J and later by R.L.Dobrushin [5] for more general deterministic 
dynamical systems. The limiting Vlasov-type equations for particle densities 
in both papers are considered in classes of integrable functions (or finite mea- 
sures in the weak form). This corresponds, actually, to the situation of finite 
volume systems or systems with zero mean density in an infinite volume. The 
Vlasov equation for the integrable functions was investigated in details by 
V.V.Kozlov [IB] . An excellent review about kinetic equations which describe 
dynamical multi-body systems was given by H.Spohn [22], [23]. Note that 
in the framework of interacting diffusions a similar problem is known as the 
McKean- Vlasov limit. 

Motivated by the study of Vlasov scaling for some classes of stochastic 
evolutions in continuum for which the use of the mentioned above approaches 
breaks down (even in the finite volumes) we developed the general approach 
to study the Vlasov-type dynamics (see [8]). It is based on a proper scaling 
of the hierarchical equations for the evolution of correlation functions and 
can be interpreted in the terms of the rescaled Markov generators. Up to 
our knowledge, at the present time only this technique may give a possibility 
to control the convergence in the Vlasov limit in the case of non-integrable 
densities which is generic for infinite volume infinite particle systems. Say- 
ing about the evolutions, the kinetic equations of which can not be studied 
by the classical techniques described in [3] and [5], we have in mind, first 
of all, spatial birth-and-death Markov processes (e.g., continuous Glauber 
dynamics, spatial ecological models) and hopping particles Markov evolu- 
tions (e.g., Kawasaki dynamics in continuum). The main difficulty to carry 
out the approach proposed by W.Braun, K.Hepp [3] and R.L.Dobrushin [5] 
for such models is absence of the proper descriptions in terms of stochastic 
evolutional equations. Another problem concerns the possible variation of 
particles number in the evolution. The important point to note also is that 
an application of the technique proposed in [8] leads to a limiting hierarchy 
which posses a chaos preservation property. 
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The aim of this paper is to study the Vlasov scaling for the individual 
based model (IBM) in spatial ecology introduced by B.Bolker and S.Pacala [TJ 
12], U.Dieckmann and R.Law jl] (BDLP model) using the scheme developed in 
[8] . A population in this model is represented by a configuration of motionless 
organisms (plants) located in an infinite habitat (an Euclidean space in our 
considerations). The unbounded habitat is taken to avoid boundary effects 
in the population evolution. 

The evolution equation for the correlation functions of the BDLP model 
was studied in details in [9] . In [TJ |2] , [I] this system was called the system 
of spatial moment equations for plant competition and, actually, this sys- 
tem itself was taking as a definition of the dynamics in the BDLP model. 
The mathematical structure of the correlation functions evolution equation 
is close to other well-known hierarchical systems in mathematical physics, 
e.g., BBGKY hierarchy for the Hamiltonian dynamics (see, e.g. [6]). As in 
all hierarchical chains of equations, we can not expect the explicit form of 
the solution, and even more, the existence problem for these equations is a 
highly delicate question. 

According to the general scheme (see [8]), we state conditions on struc- 
tural coefficients of the BDLP Markov generator, which give a weak con- 
vergence of the rescaled generator to the limiting generator of the related 
Vlasov hierarchy. Next, we may compute limiting Vlasov type equation for 
the BDLP model leaving the question about the strong convergence of the 
hierarchy solutions for a separate analysis. A control of the strong conver- 
gence of the rescaled hierarchy is, in general, a difficult technical problem. 
In particular, this problem remains be open for BBGKY hierarchy for the 
case of Hamiltonian dynamics as well as for Bogoliubov-Streltsova hierarchy 
corresponding to the gradient diffusion model. In the present paper we show 
the existence of the rescaled and limiting evolutions of correlation functions 
related to the Vlasov scaling of the BDLP model and the convergence to the 
limiting evolution. With this evolution for special class of initial conditions is 
related a non-linear equation for the density, which is called Vlasov equation 
for the considered stochastic dynamics. 

Let us mention that a version of the BDLP model for the case of finite 
populations was studied in the paper [11]. In this work the authors developed 
a probabilistic representation for the finite BDLP process and applied this 
technique to analyze a mean-field limit in the spirit of classical Dobrushin 
or McKean-Vlasov schemes. They obtained an integro-differential equation 
for the limiting deterministic process corresponding to an integrable initial 
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condition. The latter equation coincides with the Vlasov equation for the 
BDLP model derived below in our approach. 

The present paper is organized in the following way. Section 2 is devoted 
to the general settings required for the description of the model which we 
study. In Subsection 3.1 we discuss the general Vlasov scaling approach for 
spatial continuos models. Subsection 3.2 is devoted to the abstract conver- 
gence result for semigroups in Banach spaces which will be crucial to prove 
the main statements of the paper presented in Subsection 3.3. The corre- 
sponding proofs are given in Subsection 3.4. 

2 Basic fact and description of model 
2.1 General facts and notations 

Let B(M. d ) be the family of all Borel sets in M d and Bb(M. d ) denotes the system 
of all bounded sets in B(R d ). 

The space of n-point configuration is 

rf, n) = rg d :={?7CM d | \n\=n}, n e N := N U {0}, 

where \A\ denotes the cardinality of the set A. The space := 

for A e B b (R d ) is defined analogously to the space T^K As a set, is 
equivalent to the symmetrization of 

(M^ = { (xi, . . . ,x n ) e (R d ) n \ x k ^xtiik^l), 

i.e. (M. d ) n /S n , where S n is the permutation group over {l,...,n}. Hence 
one can introduce the corresponding topology and Borel cr-algebra, which 
we denote by 0(T^) and B(T^), respectively. Also one can define a mea- 
sure as an image of the product of Lebesgue measures dm(x) = dx on 
(R d ,B(R d )). 

The space of finite configurations 

To := U ^ 

neN () 

is equipped with the topology which has structure of disjoint union. There- 
fore, one can define the corresponding Borel a- algebra B(T ). 
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A set B G B(T Q ) is called bounded if there exists A G B b (R d ) and N eN 
such that B C U^Lo ■ The Lebesgue-Poisson measure X z on T is defined 

as 

V n) - 

' n! 

n=0 

Here z > is the so called activity parameter. The restriction of A 2 to Ta 
will be also denoted by X z . 
The configuration space 

r := {7 C M. d I | 7 n A| < 00, for all A G B b (R d )} 

is equipped with the vague topology. It is a Polish space (see e.g. [H]). The 
corresponding Borel a-algebra B(T) is defined as the smallest a-algebra for 
which all mappings N\ : T — >■ No, Na(j) '■= \l H A| are measurable, i.e., 

S(r)=a(iVA|AG^(R (i )). 

One can also show that T is the projective limit of the spaces {TaIae^ (R d ) 
w.r.t. the projections pa : T — > Fa, Pa(i) '■= 7a, A G £>;,(R d ). 

The Poisson measure it z on (T,B(T)) is given as the projective limit of 
the family of measures {^} A eB b (R d )j where -n z is the measure on Ta defined 
by tt^ := e~ zm ^\ z . 

We will use the following classes of functions: L° s (r ) is the set of all 
measurable functions on r which have a local support, i.e. G G L^(ro) 
if there exists A G B b {R d ) such that G fr \r A = 0; £>bs(ro) is the set of 
bounded measurable functions with bounded support, i.e. G |r \B= for 
some bounded B G B(T ). 

On T we consider the set of cylinder functions J^y^r), i.e. the set of all 
measurable functions G on (r,i3(r))) which are measurable w.r.t. B\(T) 
for some A G B b {R d ). These functions are characterized by the following 
relation: F(j) = F \ Ta ( 7a ). 

The following mapping between functions on To, e.g. L^(Tq), and func- 
tions on T, e.g. .Fcyi(r), plays the key role in our further considerations: 

KG{-y):=^2G{r,), 7 € I\ (2.1) 

where G G L^(r ), see e.g. [131 HE]- The summation in the latter 
expression is taken over all finite subconfigurations of 7, which is denoted 
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by the symbol 77 <e 7. The mapping K is linear, positivity preserving, and 
invertible, with 

K-'Fin) :=^(-l)^F(0, ^r„. (2.2) 

Let A^f m (r) be the set of all probability measures p on (T, B(T)) which 
have finite local moments of all orders, i.e. j r \'-fA\ n fi(d'j) < +00 for all 
A G B b (R d ) and n e N . A measure p on (r , B(T )) is called locally finite 
iff p{A) < 00 for all bounded sets A from B(T Q ). The set of such measures is 
denoted by A4\f(T ). 

One can define a transform K* : Aij m (T) — > .Mif(ro), which is dual to 
the ^-transform, i.e., for every p G Ai} m (T), G G i3b s (ro) we have 



KG( 1 )p(d 1 ) = / G( V )(K*p)(d v ). 
r Vr 

The measure p M := is called the correlation measure of p. 

As shown in [13] for p G .Mj m (r) and any G G L^FcbP^) the series ( 12. ID 
is p-a.s. absolutely convergent. Furthermore, KG G L 1 (r,/i) and 



G{rj) p^drj) = /(KG)( 7 )//(rf7)- (2.3) 
r </r 

A measure /i G A^J m (r) is called locally absolutely continuous w.r.t. n z iff 
Pa '■= P°P~\ 1 is absolutely continuous with respect to for all A G i3A(lR d ). 
In this case p M := K*p is absolutely continuous w.r.t X z . We denote 



The functions 



. ( R d)n R+ ( 2 4 ) 



i.(»)f T , r ^ — J ^({^ij • • • , x n}), if (xi, ■ ■ ■ ,x n ) G 

(J, otherwise 



are the correlation functions well known in statistical physics, see e.g. [20J, 
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2.2 Description of model 

We consider the evolving in time system of interacting individuals (particles) 
in the space ¥L d . The state of the system at the fixed moment of time t > 
is described by the random configuration 7t from T. Heuristically, the 
mechanism of the evolution is given by a Markov generator, which has the 
following form 

L := IT + L+ 



where 



(L-F)( 7 ) := (L-(m,H-,a-)F)(j) : = £ 



x{=7 

(L + F)( 7 ) := (L + (x + ,a-)F)( 7 ):=x^ 



m 



+ x ^ a (x-y) 



y£-y\x 



[ Y^a + {x-y)D + x F{ 1 )dx. (2.5) 



Here < a , a + G L x (]R d ) are arbitrary, even functions such that 

a~(x)dx = / a + (x)dx = 1 



(in other words, a", a + are probability densities) and m, x~, x + > are 
some positive constants. 

The pre-generator L describes the Bolker-Dieckmann-Law-Pacala BDLP 
model, which was introduced in [HEJH]. During the corresponding stochastic 
evolution the birth of individuals occurs independently and the death is ruled 
not only by the global regulation (mortality) but also by the local regulation 
with the kernel x~a~. This regulation may be described as a competition 
(e.g., for resources) between individuals in the population. 

The evolution of the one dimensional distribution for such systems can 
be expressed in terms of their characteristics, e.g. the correlation functions 
(see f)2.4p ). The dynamics of correlation functions for the BDLP model was 
studied in [9]. The main result of this paper informally says the following: 

If the mortality m and the competition kernel xTa~ are large enough, 
then the dynamics of correlation functions, associated with the pre-generator 
(12. 5ft . exists and preserves (sub-)Poissonian bound. 

For the readers convenience we repeat below the relevant material from 
j9] without proofs. 
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Let L ± := K^L^K be the K-image of L ± , which can be initially defined 
on functions from .BbstTo)- For arbitrary and fixed C > we consider the 
operator L := L + + L~ in the functional space 

c c = L 1 (r ,cWd\ (77)). 

Below, symbol stands for the norm of this space. 

For any u>0we define H{u) to be the set of all densely defined closed 
operators T on C C) the resolvent set p(T) of which contains the sector 

Sect := {CeC J |argC| < \+^), 

and for any e > 

IKT-am^H, |argC|<^ + ^-e, 

where M £ does not depend on (. 

The first non-trivial result, which is based on the perturbation theory, 
says that the operator L with the domain 

D(L) :=[GeCc\ \-\G(-)eC c , E a ~ (•)£(•) G C c ) 

is a generator of a holomorphic C -semigroup U t on £ c . 

To construct the corresponding evolution of correlation functions we note 
that the dual space {Cc)' = (L 1 (r , rfAc)) = -^(To? d\c), where d\c '■= 
C^dX. The space {Cc)' is isometrically isomorphic to the Banach space 

k c ■= [k : r -)• r /c(-)c _ ' ' g L°°(r , A) I 

with the norm 

\\k\\ Kc := ||C7-Hfc(.)||L-cro,A), 
where the isomorphism is provided by the isometry Rc 

{C c )' 3 k ^ R c k := k{-)C 11 G/C c . (2.6) 

In fact, we have duality between Banach spaces Cc and JCc given by the 
following expression 

{{G, k)) := / G-kdX, GeC c , keJCc (2.7) 
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with 

\((G,k))\<\\G\\ c -\\k\\ Kc . (2.8) 
It is clear that for any k G /Cc* 

\k{rj)\<\\k\\ Kc C ]r >\ for A -a.a. 7] G r . (2.9) 

Let L' be the adjoint operator to L in {Cc) 1 with domain D(L'). Its image 
in K-c under the isometry R c we denote by L* = R c L'R c -i- It is evident 
that the domain of L* will be D(L*) = RcD(L'), correspondingly. Then, for 
any G G C c , k G D(L*) 

G-L*kd\= I G ■ R c L'Rc-ikdX — I G-L'R c -ikd\ c 

To J To J To 



LG-R c -ikd\ c = / LG-kdX, 

To JTo 

therefore, L* is the dual operator to L w.r.t. the duality (12. 7ft . By [10J, we 
have the precise form of L*\ 

(L*k)(rj) = - (m\r)\ + x~E a ~ (77)) k(rj) (2.10) 

+ x+ ^2 ^2 a+ ( x - y)Hv \ x) 

x£rj y£r)\x 



+ x-/Va-( I -#(,U l)t 



Now we consider the adjoint semigroup T'(t) on (Cc)' and its image T*(t) in 
/Cc*. The latter one describes the evolution of correlation functions. Trans- 
ferring the general results about adjoint semigroups (see, e.g., [7J) onto 
semigroup T*(t) we deduce that it will be weak*-continuous and weak*- 
differentiable at 0. Moreover, L* will be the weak*-generator of T*(t). Here 
and subsequently we mean "weak*-properties" w.r.t. the duality (12. 7p . 
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3 Vlasov scaling 

3.1 Description of Vlasov scaling 

We begin with a general idea of the Vlasov-type scaling. It is of interest to 
construct some scaling L e , e > of the generator L, such that the following 
scheme works. 

Suppose that we know the proper scaling of L and we are able to prove 
the existence of the semigroup T e (t) with the generator L e := K~ x L e K in the 
space Cc for some C > 0. Let us consider the Cauchy problem correspond- 
ing to the adjoint operator L* and take an initial function with the strong 
singularity in e. Namely, 

$\r))~e- M ro(r}), e -> 0, rj E Tq, 

where the function r$ is independent of e. The solution to this problem is 
described by the dual semigroup T*(t). The scaling L i-> L e has to be chosen 
in such a way that T*(t) preserves the order of the singularity: 

Another very important requirement for the proper scaling concerns the dy- 
namics ro i — y r t . It should preserve Lebesgue-Poisson exponents: if ro(r/) = 
e\(po, rj) then r t (r]) = e\(pt, v) an d there exists explicit (nonlinear, in general) 
differential equation for p t 

d 

— pt(x) = v(pt(x)), (3.1) 

which will be called the Vlasov-type equation. 

Now let us explain the main technical steps to realize Vlasov-type scaling. 
Let us consider for any e > the following mapping (cf. fl2.6j> ) on functions 
on r 

(R £ r)(v) :=e^r(ri). (3.2) 

This mapping is "self-dual" w.r.t. the duality (12. 7p . moreover, Rj 1 = R e -i. 
Then we have ~ R E -ir , and we need r t ~ R £ T*(t)k^ ~ R £ T*(t)R £ -ir . 
Therefore, we have to show that for any t > the operator family R e T*(t)R e -i , 
£ > has limiting (in a proper sense) operator U (t) and 



U{t)e x (p ) = ex(pt). 



(3.3) 
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But, informally, T*(t) = exp {tL*} and R £ T*(t)R £ -i = exp {tR £ L*R £ -\} . Let 
us consider the "renormalized" operator 

Ll ren :=R £ t £ R £ -,. (3.4) 

In fact, we need that there exists an operator V* (called Vlasov operator) such 
that exp {tR £ L* £ R £ -i} -»■ exp {tV*} =: U(t) for which <g£S§ holds. Hence, 
heuristic way to produce the scaling L i— > L £ is to demand that 

Jim (J^e x (p t ,ri) ~ L* £)Tea e x {p t , =0, 77 G T , 

if p t satisfies (13. ip . The point-wise limit of L* ren will be natural candidate for 
V*. Having chosen the proper scaling we proceed to the following technical 
steps which give the rigorous meaning to the idea introduced above. Note 
that definition (13 .4p implies L £ ^ cn = R £ -\L £ R £ . We prove that "renormal- 
ized" operator L £>ren is a generator of a contraction semigroup T £ ^ cn (t) on 
Cc- Next we show that this semigroup converges strongly to some semigroup 
Ty(t) with the generator V. This limiting semigroup leads us directly to the 
solution for the Vlasov-type equation. Below we show how to realize this 
scheme in details. 



3.2 Approximation in Banach space 

In this subsection we study general question about the strong convergence 
of semigroups in Banach spaces. The obtained results will be crucial in the 
realization of the Vlasov-type scaling for the BDLP model. 

Let {Ut, t > 0} , e > be a family of semigroups on a Banach space 
E. We set (L £ , D(L £ )) to be the generator of {Uf , t > 0} for each e > 
0. Our purpose now is to describe the strong convergence of semigroups 
{Uf, t > 0} , e > in terms of the corresponding generators as e tends to 
0. According to the classical result (see e.g. [12]). it is enough to show that 
there exists (3 > and A : Re A > (3 such that 

(L £ - Al)" 1 (L - Al)" 1 , e^0 } (3.5) 

where 1 is the identical operator. In this subsection we show how to prove 
(13. 5 P under the following assumptions on the family (L £ , D(L £ )), e > 0: 
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Assumptions (A): 

1. For any e > 0, the operator (L e , D(L £ )) admits representation 

L e = Ai(e) + A 2 (e), 
where A\{e) is a closed operator and D(Ai(e)) = D(A2(e)) := D(L £ ). 

2. There exists > and A: Re A > /3 such that 

(a) A belongs to the resolvent set of Ai(e) for any e > and 
(A x (e) - Xty 1 (Ai(0) - Al)" 1 ,e ->• 0, 

(b) 

sup - Al) _1 || < ||(Ai(0) - AH)" 1 !! , 

£>0 

(c) for any e > 

||A 2 (e) (Ax^) - AH)" 1 !! < 1, 

(d) (A 2 (e) (Ai(e) — Al) -1 + l) 1 converges strongly to the operator 
(A 2 (0) (Ai(0) - Al) -1 + as e -> 0. 

The strong convergence result for the family {[/ t e , t > 0} , e > is estab- 
lished by our next theorem 

Theorem 3.1. Let (L e , D(L e )), e > be the family of generators corre- 
sponding to the family of C -semigroups {Uf, t > 0} , e > 0. Then, Uf con- 
verges strongly to Z7 t ° as £ — > uniformly on each finite interval of time, 
provided assumptions (A) are satisfied. 

Proof. The proof is completed by showing ( 13. 5p . For any e > and A from 
the resolvent set of A\{e) we have 

Ran ((Ai(£) - Al)- 1 ) = D {A 1 {e)) = D {A 2 (e)) = D{L £ ). 

Hence, 

L e -Xl = A x {e)+A 2 (e) - Al 
= (A a (e) (^(e) - Al)" 1 + l) {A x {e) - Al) . (3.6) 



Operator approach to Vlasov scaling for models of spatial ecology 



13 



Combining (13. 6 p with the assumption 2(c) of (A) we get the following rep- 
resentations for the resolvent 



(L £ - Xiy 1 = (Ax(e) + A 2 (e) - XI)' 1 
(A 1 (e) - Xiy 1 (A 2 (e) (A 1 (e) - Xiy 1 + l)' 1 . 



(3.7) 



From this formula, triangle inequality and assumptions 2(a), 2(b) and 2(d) 
of (A) we conclude the assertion of the theorem. □ 



3.3 Main results 

We check at once that the proper scaling for the BDLP pre-generator is the 
following one 



+ X 



ni 



+ ex a {x — y) 



D-F^) 



(3.8) 



f Y^a + (x~y)DtF( 1 )dx, 



e > 0. 



Next we consider the formal i^-image of L £ and the corresponding renor- 
malized operator on £? bs (r ): 

L e G := K~ 1 L £ KG; L £jren G := R e -iL e R E G. 

In the proposition below we calculate the precise form of the operator L £j ren 
for the BDLP model. 

Proposition 3.2. For any e > and any G e B^ s (Tq) 

L £ ,renG =A X G + A 2 G + e (B 1 G + B 2 G) , 
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where 



(A X G) (77) = — m\rj\G (77) , 

(A 2 G) (77) = - x~ ^ £ a- (x-y)G(iAx) 



+ x+E/ , « + (»-»)Gta\«u*)«fa, 



yet] 

(B 1 G) {r,)=-H-E a - ( V )G(rj), 

(B 2 G) (77) = x + V / a + (x - y) G (77 U x) cfx. 
ycv jRd 

Proof. According to the definition, we have Z/ £jren = R E -iL e R e , where 
L e = L~ (m, ek~cl~) + e~ 1 L + (ex + a + ) . 

As a result, 

(L £ G) (77) = {A\G) (77) + e{BxG) (77) + (B 2 G) (77) 
- a_ ( s - Z/) G (^/ \ ^) 

x£r/ y£rj\x 



+ „+(x- w )G(,\ w U,)«fc. 



and hence 

(L £ ,renG) (77) = (A X G) (77) + (A 2 G) (77) + £((£?! + S 2 ) G) (77) , 

which completes the proof. □ 

Remark 3.3. It is easily seen that the operator V := A\ + A 2 will be the 
point- wise limit of L £)Ten as e tends to 0. Therefore, the adjoint operator to 
V w.r.t. to the duality (12. 7\i (if it exists) can be considered as a candidate 
for the Vlasov operator in our model. 

Below we give the rigorous meaning to the operator L £iXen . Let us define 
the set 

D\ := [G g C c \ E a ~ (■) G (•) G C c , H G (•) G £ c } 
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Proposition 3.4. For any e, m, x , C > the operator 

Axis) := A x + eB 1 (3.9) 

with the domain D x is a generator of a contraction Co-semigroup on Cc- 
Moreover, A x {e) e H (u) for all u G (0; f ) . 

Proof. See the proof of Proposition 4.2 in [9]. □ 

Remark 3.5. It is a simple matter to check that Proposition 13.41 holds also in 
the case e = 0, provided the domain of the operator Ai(0) := A x is changed 
to 

D := {G e C c | \-\GeC c }DD 1 . 
The next task is to show that for any e > the operator 

A 2 {e) := L £ ,ren - A^e) = A 2 + eB 2 (3.10) 

with the domain D x as well as the operator A 2 (0) := A 2 with the domain 
D are relatively bounded w.r.t. the operator (Ai(e), D x ) and (A x , D ), 
correspondingly. This is demonstrated in Propositions 13.61 and 13. 7\ which 
can be proved similarly to Lemmas 4.4 and 4.5 in [9]. 

Proposition 3.6. For any 5 > and any x~ , x + ,m,C > such that 

h-C x+ . 

+ — < 8 

m m 

the following estimate holds 

\\A 2 G\\ C ^SWA&Wc, G G Do. 
Moreover, for all e > 

\\A 2 G\\ c <5\\A 1 {e)G\\ c , G e D x . 
Now, the operator (A 2 , D ) is well-defined on Cc- 

Proposition 3.7. For any e, 5 > and any x~, x + , m, C > such that 

ex + E a+ {n) < 5C (ex E a ~ (77) + m \r)\\ , rj ^ 
i/ie following estimate holds 

\\eB 2 G\\ c < a ||A 1 (e)G|| c , G e D x 

with a < 5. 
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Remark 3.8. Proposition 13.71 enables us to take D{B<i) = D\. As a result, 
Remark 13.51 shows that the domain of the operator Az(e) will be Dq fl D\ = 

We are now in a position to show that the operator (L e ^ ien , Di) generates 
semigroup on Cc- To this end we use the classical result about the pertur- 
bation of holomorphic semigroups (see, e.g. [TJ]). For the convenience of the 
reader we formulate below the main statement without proof: 

For any T G T-L{u), u G (0; |) and for any e > there exist positive 
constants a, 6 such that if the operator A satisfies 

\\Au\\ < a\\Tu\\ +b\\u\\, u G D(T) c D(A), 

with a < 5, b < 5, then T + A is a generator of a holomorphic semigroup. In 
particular, if b = 0, then T + A G H(oj — e). 

Theorem 3.9. Let the functions a", a + and the constants m, x~, x + , C > 

satisfy 

m> A{k~C + x + ) , (3.11) 
Ck~oT (x) > 4x+a+ (x), x G R d . (3.12) 

Then, for any e > the operator (L £yTen , Di) is a generator of a holomorphic 
semigroup U t , e , t > on Cc- 

Proof. Let e > be arbitrary and fixed. By definition, 

L 6 ,r 8 n = ^i(e) + ^ 2 (e). 

The direct application of the theorem about perturbation of holomorphic 
semigroups (see the formulation above the assertion of Theorem 13.91) to T = 
Ai(e) and A = Ai{e) gives now the desired claim. It is important to note 
that Proposition 13 .41 enables us to consider 5 equal to | in the formulation of 
the classical theorem introduced above. The appearance of the multiplicand 
4 on the left-hand side of the both assumptions in assertion of Theorem 13.91 
is motivated exactly by the latter fact. □ 

Theorem 3.10. Assume that the constants m, x~ , x + ,C > satisfy 

m > 2 (x^C + x + ) . 

Then, the operator V = A\ + A 2 with the domain D is a generator of a 
holomorphic semigroup UY, t > on C c . 
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Proof. We use the same classical result as for Theorem 13.91 in the case: A\ 
is a generator of holomorphic semigroup, A2 is relatively bounded w.r.t. A\ 
with the boundary less then |. □ 

Now we may repeat the same considerations as at the end of Section [2j 
Namely, transferring the general results about adjoint semigroups (see, e.g., 
[7]) onto semigroup {UjY m we deduce that it will be weak*-continuous 
and weak*-differentiable at 0. Moreover, V* will be the weak*-generator of 
T*(t). This means, in particular, that for any G G D(V) C Cc, k G D(V*) C 

j f ((G, (Uryk)) = ((G, V*(U?)*k)) . (3.13) 
The explicit form of V* follows from (I2.10p . namely, for any k G -0(1^*) 



V*k(i]) = — m\n\k(rj) 



x 



+ I a+ (x -y)k(rj\xUy)dy. (3.14) 




As a result, we have that for any k G D(V*) the function k t = (U^)*k 
provides a weak* solution of the following Cauchy problem 



(3.15) 



In the next theorem we show that the limiting Vlasov dynamics has chaos 
preservation property, i.e. preserves the Lebesgue-Poisson exponents. 

Theorem 3.11. Let conditions of Theorem \3.9i be satisfied and, additionally, 
C > 16 g_ 1 ■ Let p > be a measurable nonnegative function on M d such that 
esssup xgR d po(x) < C . Then the Cauchy problem (13.151) with k$ = e\(po) has 
a weak* solution k t = e\(pt) G JCc, where p% is a unique nonnegative solution 
to the Cauchy problem 

d 

—p t (x) = x + {a + * p t )(x) - x p t (x)(a * p t )(x) - mp t (x), 



and esssup^gRd pt(x) < C , t > 0. 



Operator approach to Vlasov scaling for models of spatial ecology 



18 



Proof. First of all, if (13. 161) has a solution pt{x) > then 
d 

— p t (x) < x + (a + * p t )(x) - mp t (x) 

and, therefore, pt{x) < r t (x) where r t (x) is a solution of the Cauchy problem 
d 

— r t (x) = x+(a + * r t )(x) - mr t (x), 

, r 4=o( x ) = Po( x ) > °> 
for G R . Hence, 

r t (x)=e-( m - + V + ^+p (x), 

where 



(L a+ f)(x):= / a + (x-y)[/(y)-/(x)]dy. 

Since for / G L°°(R o! ) we have |(L 0+ /)(o;)| < 2||/|| Xoo(Kd) then, by (J3TU), 

r t (x) < (7 e -(^-^ + )* e 2^ + * < c, 

that yields < p t (x) < C. 

To prove the existence and uniqueness of the solution of ( 13.161) let us 
fix some T > and define the Banach space X T = C([0; T], L°°(R d )) of all 
continuous functions on [0; T] with values in L°°(R d ); the norm on Xt is given 
by ||m||t := max ||wt||nx>fRdv We denote by X$ the cone of all nonnegative 

te[0;T] K ' 

functions from X?. 

Let $ be a mapping which assign to any v G Xt the solution u t of the 
linear Cauchy problem 

d 

— u t (x) = x+(a+ *v t )(x) - x u t (x)(a * v t )(x) - mu t (x), (3 17) 

for G R . Therefore, 

($f ) t (x) =exp |— J (m + x~(a~ * f s )(x))<is| po(x) (3.18) 

+ / exp < — / (m + x~(a~ * v T )(x))dr > x + (a + * v s )(x)ds. 
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We have that v £ implies $w > as well as the estimate 



x' 



(®u)t(x) < p (x) + x + \\v\\ T / e- {t - s)m ds <C+ — i|v|| T , 

Jo ™ 

where we use the trivial inequality 

\\f*g\\ L ~m<\\f\\ L HR*)\\9h°°m> /e^W- geL°°(R d ). (3.19) 
Therefore, $w £ X^. For simplicity of notations we denote for v £ 

(Bv)(t, x) = m + x~{aT * v t )(x) > m > 0. 
Then, for any v, w £ X£ 

\($v)t(x) - ($w)t(x)\ 



< 



exp < — / (Bv)(s, x)ds > — exp < — / (Bw)(s,x)ds 



Po{x) 



+ 



exp 



— exp 



(Bv)(t, x)dr > x + (a + * v s )(x) 



(Bw)(t, x)dr > x + (a + * w s )(x) 



ds. 



We have 



<e 



<e 



exp <^ — / (Bv)(s, x)dsj — exp |— J (Bw)(s,x)ds 

x~(a~ * v s )(x)ds\ — exp |— J x~(a~ * w s )(x)ds 

t ,t 
x (a *v s )(x)ds— / x (a *w s )(x)ds 



-rat 



-rat 



exp 



o 



<e mt x \\v-w\\ T -t< — \\v — w\\ T , 

em 

where we used (I3.19P and obvious inequalities \e~ a — e~ b \ < \a — b\ for a, b > 0; 
e~ x x < e -1 for x > 0. 

Next, using another simple estimates for any a, b,p, q > 

\pe- a - qe~ b \ < e' a \p - q\ + qe- b \e~ {a - b) - 1| < e - a \p -q\ + qe~ b \a - b\, 



x 
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(Bv)(t, x)dr > x + (a + * v 8 )(x) 



exp 



— exp < — / (Bw)(t, x)dr > x + (a + * w s )(x) 



ds 



<x + J exp |— J (Bv)(t, x)dr^ \a + * (v s — w s )\(x)ds 
+ / exp < — / (Bw)(t, x)dr > (n + a + * w s )(x) 



x 



t ft 
(Bv)(r,x)dT— / (Bw)(T,x)dr 



- m(t - s) ds 



ds 



<x \\v — w\\t I e 
'o 

+ J exp |— J K~{aT * w T )(x)dr^ (x + a + * w s )(x) 
x e~ m{t ~ s ^ I x~(cT * \v T - w T \)(x)drds 



and, using f 1 3 . X 2 j) and the inequalities above, one can continue 



X 1 

< — \\v — w\\t H lit; — w\\t 

m 4 em 

-t ( rt 



x / exp 



o 



X 



(a *w T )(x)dr>>c (a * w s )(x)ds 



x+„ , Cx", 

: — u — w hp + \\v — W t 

m 4 em 

«t 



< 







d 


X 


/' 

/o 


— exp 

as 


x+ 


+ 


Cx~\ 






m 




4 em / 



/ ^-expi— I x (a * w T )(x)dr\ ds 
Jo ds { J s J 



\V -10 y. 



Therefore, for u , u> G X~ 
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if, e.g., 1 + j < 4Ce, that means C > 16 ^_ 1 . 

As a result, by (13. lip . $ is a contraction mapping on the cone X£ . Taking, 
as usual, w (n) = $V 0) , n > 1 for «(°) G X+ we obtain that {w (n) } C X+ 
is a fundamental sequence in Xt which has, therefore, a unique limit point 
v G Xt- Since X^~ is a closed cone we have that v G XyT . Then, identically to 
the classical Banach fixed point theorem, v will be a fixed point of $ on Xt 
and a unique fixed point on X£. Then, this v is the nonnegative solution of 
(I3.16P on the interval [0;T]. By the note above, v t (x) < C. Changing initial 
value in ( 13.161) onto pt\ t=T {x) = vt(x) we may extend all our considerations 
on the time- interval [T; 2T] with the same estimate v t (x) < C; and so on. 
As a a result, (13. 16ft has a unique global bounded non- negative solution pt(x) 
on R+. 

Consider now 

hiv) = e\{Pt,v) ^ K, c , 

Using (13.161) and (I3.14p . we immediately conclude that k t {rj) = e\(p t ,rj) is a 
solution to (I3.15p . □ 

The main result of the paper is formulated in the next theorem. Its proof 
will be given in Subsection 3.4. 

Theorem 3.12. Under conditions of Theorem \ 3.9i the semigroup U t , £ con- 
verges strongly to the semigroup as e — > uniformly on any finite inter- 
vals of time. 



3.4 Proofs 

According to Theorem l3.ll the statement of Theorem 13 . 1 21 will be proved once 
we verify Assumptions (A) for the operators (Ai(e), Di), (^(e), Z?i), e > 0, 
defined in the previous subsection. Note, that Ai(0) = A\ and A 2 (0) = A 2 
are defined on the domain Dq. 

In the following proposition we verify Assumption 2(a) of (A). 

Proposition 3.13. Let A > then 

{Axis) - AH)" 1 (A x - \\y l ,e 0. 
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Proof. For any G G £ 



c 



(Ai(e) - Aip G — (Ai — Al) -1 G\ 

1 



c 



r 



\—m \rj\ — e>c~E a ~ (77) — A — m \ 
\G(r,)\F e (r,)CMd\(r,), 



C M dX (77) 



where 



exT E a {rj) 



(m\rj\ + ex~E a (rj) + A) (m \rj\ + A) 



Since < F £ (77) < 1/A and F £ (rj) — > as e — > for any 77 6 r , we get the 
desired statement. □ 

Next we check Assumption 2(b) of (A). 

Proposition 3.14. Let A > be arbitrary and fixed. Then 

sup ||(Ai(e) - All) -1 1| < ||(Ai - Al)~ 



x — II 



£>0 



Proof. For any G G Cc an d an Y £ > 
||(Ai(e)-Al) -1 G" 
\G(v)\ 



\c 



< 



m\rj\+ex~E a (rj) + A 
1 



C |r?l dA (77) 



\\G\\ 



C^dX{r]) = ||(Ai - A1) _1 G| 



v — II 



c ■ 



< ||(Ai - Al)~ 

This finishes the proof. 

Assumption 2(c) of (A) is proved in the next Proposition. 



□ 



Proposition 3.15. Let conditions of Theorem \ 3.9\ be satisfied. Then, for 
any A > 

sup||A 2 (e)(A 1 (e)-Al)- 1 || <\ (3.20) 



£>0 
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Proof. First we prove assertion for e — 0. Since D (Ax) = D (A 2 ) = Dq and 
Ran ((Ax - Al)" 1 ) = D (Ax), the operator A 2 (A 1 - AH) -1 is well defined. 
Next, inequality (13. lip and Proposition 13.61 yields 

||A 2 (A a - Al) -1 !! < |. (3.21) 

Indeed, 

\\A 2 G\\ c <a\\A 1 G\\ c <a\\(A 1 -Xl)G\\ c 
with a < |. Therefore, 

||A 2 (A 1 -A1)- 1 G|| C < 

and (I3.2ip is proved. 

Now, let £ > be arbitrary and fixed. The main arguments we use to 
show 

\\A 2 (e) {A 1 {e)-Xl)- 1 \\<^ 

are the following: 

1) D(A x (e)) = D x C D = D(A 2 ). Hence, A 2 (A x (e) - Al) -1 is well- 
defined on Cc- Moreover, Proposition 13.61 implies 

||A 2 (A 1 (e)-Al)" 1 || < i e>0. 

2) D (B 2 ) = D (A x (e)) = £>i and for any e > 

||eB 2 (Ai(e)-Al)- 1 || < 

which follows from Proposition 13.71 

3) Since A 2 (e) := A 2 + e_B 2 , we have 

||A 2 (e)(A 1 (e)-Al)- 1 || < ~ (3.22) 

The latter concludes the proof. □ 
We set 

Q £ = (A 2 {e) (Axis) - Al)" 1 + l)" 1 , Q = (A 2 (Ax - XI)' 1 + I)" 1 . 
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In order to verify Assumption 2(d) of (A) we have to show that Q £ — > Q 
as e — > 0. 

Suppose that we can show that 

A 2 (A 1 (e) - Xiy 1 A 2 (A, - Xiy 1 , s 0. 
eB 2 {A l {e) -Xiy 1 -^0, e^O. "~ 0> 

Then, 

C £ :=A 2 (e) (A^-Xiy 1 

=A 2 {A l + eB x - A!)" 1 + sB 2 {A x + eB x - Xl)' 1 A 2 {A x - Xiy 1 



To check 



we proceed as follows: 

i^ + iy'-Q 



q £ = (c £ + ty 1 g (3.24) 



= (C £ + l)" 1 - (A 2 (A, - xiy 1 + i) _1 

= (C 6 + I)" 1 (A 2 (Ax - Al)" 1 + 1 - C £ - 1) (A 2 (Ax - Al)" 1 + l)~ 1 
= (C e + I)" 1 (A 2 (Ax - Al)" 1 - C e ) (A 2 (Ax - Al)" 1 + l) _ 1 . 

Assuming (I3.23f) it is obvious now that convergence (13. 24f) is equivalent to 



sup II (C £ + 1) 1 1| < oo, 

£>0 



which is clear from 



Ce + ir 1 !! < x.j^n and lic-ii < |- 

The last bound we conclude from (I3.22p . As a result we shall have established 
Theorem E12 if we show f l3T23|) . 

Lemma 3A6. A 2 {A x {e) - Al) -1 A 2 (Ax - Al) -1 , as £ 0. 
Proof of Lemma \3.1b\ Proposition 13.61 and 

D (Ax( £ )) =J}xCD (Ax) = £> (A 2 ) = D 
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leads to the following formula 

A 2 (A l (e) - Xiy 1 = A 2 (A, - Al)" 1 (A, - Al) (A^e) - AH)" 1 
Now, we are left with the task to show that 

(A l - Al) (Ax{e) - Al)" 1 -A 1, as e -> 0. 
But, for any G G Lq 

|| ((Ai - Al) (Ai(e) - Al)" 1 - 1) G 

1 



ic 



r 



m 


r/ 


+ A 


m 


»7 


+ ex" 


£ a ~ (r?) + A 



£x-,E a (77) 



G(n)\C^dX (rf) 
G{r])\C lril dX{n) -> 0, as £ 



/ r m |?y| + ex~E a (77) + A 

due to the Lebesgue's dominated convergence theorem 
Lemma 3.17. e£ 2 (Ai(e) - Al) 



-1 s 



0, as £ — > 0. 



Proof of Lemma \3. 1 7[ Since ||S 2 G|| C < \ \B\G\ C , we have to show that 
||eSi(Ai(e) - X\y x G\ c ^ 0, as e ->■ 0. 

But, 

||eSi (Ai(e) - X\y x G\ c 

= 1 gx W _ |G(r7)|C |T?l dA(r7) ->• 0, £ -)> 0. □ 

y^mlr/l+ex-^-^ + A 1 W;| KU 

The last two lemmas conclude the proof of the main Theorem. □ 



Remark 3.18. Under assumptions of Proposition 13.151 we get the following 
representation for the resolvents of V and L £ ren 

(V - Xt^' 1 = (A 1 + A 2 - Xiy 1 = - Al)" 1 (A 2 (Ai - Al)" 1 + l)~ 1 , 



I . \ II ] =(Ai(e)+A 2 (e)-Al) 1 (3.25) 
= (Ai(e) - Al)" 1 (A 2 (e) (A^e) - Al)" 1 + l)" 1 , A > 0. 
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